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ABSTRACT:  An  in  silico  screen  of  the  NIH  Molecular  Library  Small 
Molecule  Repository  (MLSMR)  of  ~350000  compounds  and 
confirmatory  bioassays  led  to  identification  of  chaetochromin  A  (l) 
as  an  inhibitor  of  botulinum  neurotoxin  serotype  A  (BoNT  A). 
Subsequent  acquisition  and  testing  of  analogues  of  1  uncovered  two 
compounds,  talaroderxines  A  (2)  and  B  (3),  with  improved  activity. 
These  are  the  first  fungal  metabohtes  reported  to  exhibit  BoNT/A 
inhibitory  activity. 
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Botulinum  neurotoxins  (BoNTs)  are  the  etiologic  agents  of 
the  neuroparalytic  disorder  botuhsm.  BoNTs  are  among 
the  most  potent  toxins  known  and  are  classified  by  the  Centers 
for  Disease  Control  and  Prevention  as  one  of  the  six  highest 
risk  threat  agents  for  bioterrorism.  The  limitations  of  current 
therapies  against  botulism  have  prompted  research  for 
alternative  therapeutics  that  could  be  used  to  treat  pre-  or 
postexposure  BoNT  intoxication,  including  studies  to  discover 
and  develop  small  molecular  weight  inhibitors  of  BoNT. 

After  demonstrating  that  a  strategy  of  in  silico  screening, 
followed  by  tiered  bioassay  evaluation  of  screening  hits,  could 
successfully  lead  to  confirmed  inhibitors  of  botulinum  neuro¬ 
toxin  serotype  A  (BoNT/A),'  we  conducted  a  new  round  of  in 
silico  screening,  this  time  using  the  NIH  Molecular  Library 
Small  Molecule  Repository  (MLSMR)  collection^  of  ~350000 
compounds,  for  structures  with  likely  binding  affinity  for  the 
active  site  of  BoNT/A  Among  the  first  hits  identified  and 
confirmed  active  from  that  screen  was  chaetochromin  A  (l),  a 
bis-naphtho-y-pyrone  produced  by  many  species  of  Chaeto- 
mium^  and  more  recently  found  in  Fusarium  spp.’^  Herein,  we 
report  details  of  our  in  silico  screen  and  the  evaluation  of 
chaetochromin  A  and  several  related  compounds  (2—6). 
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The  in  silico  screen  (Figure  SI  in  the  Supporting 
Information)  of  the  MLSMR  was  performed  using  Docking- 
based  Virtual  Screening  Pipefine  (DOVIS).^  DOVIS  was  used 
to  systematically  dock  the  MLSMR  compounds  to  the  catalytic 
binding  site  of  the  BoNTA  protein  target  (PDB  code  3C8A).^ 
The  compounds  were  sorted  by  AutoDock  4.0  score,^  and  the 
top  20000  compounds  were  chosen  and  minimized  within  a 
fixed  BoNT/A  catalytic  site.  After  minimization,  the  com¬ 
pounds  were  rescored  with  three  different  scoring  functions, 
AutoDock  4,^  LigScore  2,®  and  X-Score.^  Three  rescored  fists 
were  generated,  with  compounds  sorted  by  corresponding 
docking  score  within  each  fist.  The  top  500  compounds  from 
each  fist  were  pooled  for  further  analysis. 

Two  fists  of  compounds  were  generated  at  this  point.  The 
first,  the  consensus  fist,  was  populated  with  300  compounds 
that  appeared  in  aU  three  rescored  hit  fists.  The  second,  the 
MET-TOX  fist,  consisted  of  570  compounds  that  passed  an 
imposed  metabolism-toxicity  (MET-TOX)  filter,  which  con- 
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tained  mutagenic/®  hepatotoxic/^’^^  and  CYP-2D6^^’^* 
components.  The  consensus  and  MET-TOX  lists  were 
compared,  and  all  unique  compounds  were  retained  in  a 
combined  list.  The  surviving  compounds  were  visually 
inspected  in  the  BoNT/A  binding  site.  Using  prior  experience 
with  this  target,  we  removed  compounds  with  unrealistic 
interactions  or  geometries,  in  particular,  compounds  incapable 
of  coordination  with  the  active  site  Zn^^. 

Initially,  a  total  of  100  chemically  diverse  compounds  were 
selected  from  the  combined  list  for  testing  (Figure  S2  in  the 
Supporting  Information)  for  inhibition  of  the  protease  activity 
of  BoNT/A  LC  using  an  HPLC-based  assay.^  Chaetochromin 
^3,4,19,20  (^]y[L5000563176,  l)  was  found  to  be  the  most  active 
inhibitor  in  the  group  tested.  A  sample  of  chaetochromin  A 
obtained  from  the  NCI  Repository  (NSC  345647)  for 
confirmatory  testing  exhibited  only  ~50%  inhibition  at  100 
^M,  but  subsequent  HPLC-ELSD-MS  analysis  of  this  sample 
showed  that  it  was  <85%  pure. 

Samples  of  pure  (>95%  by  NMR)  chaetochromin  A  (l)  and 
several  related  compounds — talaroderxines  A  and  B  (2  and  3), 
secalonic  acid  A  (4),  cephalochromin  (5),  and  skyrin  (6) — 
were  obtained  from  the  fungal  compound  collection  at  the 
University  of  Iowa  for  comparative  testing.  Chaetochromin  A 
(l)  was  isolated  from  solid-substrate  fermentation  cultures  of 
Chaetomium  arcuatum  (syn.  Chaetomium  virescens)  (NRRL 
25243  =  IMI  86456)  isolated  from  a  soil  sample  collected  in 
Lucknow,  India.^^  Talaroderxines  A  and  B  (2  and  3)  were 
obtained  from  liquid  cultures  of  a  coprophilous  Delitschia  sp. 
(JS  300)  isolated  from  a  sample  of  kangaroo  dung  collected  in 
Australia.  The  talaroderxines  were  originally  described  as 
metabohtes  of  Talaromyces  derxii  cultured  from  cultivated  soil 
in  Kurashiki,  Japan,^^  and,  to  our  knowledge,  have  not  been 
previously  reported  from  any  other  source.  Secalonic  acid  A 

(4)23,24 

was  obtained  from  a  culture  of  Setophoma  terrestris 
(NRRL  25008)  obtained  as  a  colonist  of  an  Aspergillus  flavus 
sclerotium  that  was  buried  in  sandy  field  soil  near  Kilbourne, 
IL.  Cephalochromin  (5)i®>20,25,26  obtained  from  a  culture 
of  Cosmospora  sp.  (syn.  Acremonium  butyri)  (NRRL  28291) 
that  was  isolated  as  a  colonist  of  the  stromata  of  Hypoxylon  sp. 
found  on  a  hardwood  log  in  Peoria,  IL.  Skyrin  (6)^^”^^  was 
isolated  from  a  culture  of  Geomyces  pannorum  (NRRL  22978) 
obtained  as  a  colonist  of  an  Aspergillus  flavus  sclerotium  that 
was  buried  in  sandy  field  soil  near  Kilbourne,  IL.  Compounds 
of  this  general  type  usually  display  axial  chirality  (i.e., 
atropisomerism)  because  of  the  high  energy  barrier  to  rotation 
at  the  bond  linking  the  individual  aromatic  systems  in  each 
structure.  The  axial  configurations  of  1—3,  5,  and  6  have  been 
assigned  and  are  reflected  in  the  structures  shown  (Chart  l). 

Table  1  summarizes  the  results  for  both  our  primary  HPLC- 
based  protease  assay  and  the  tissue-based  mouse  phrenic  nerve 
hemidiaphragm  assay  (MPNHDA).  Not  only  was  the  activity  of 
chaetochromin  A  (l),  suggested  by  the  in  silico  screen, 
confirmed,  but  talaroderxines  A  (2)  and  B  (3)  proved  to  be 
more  active  in  the  protease  assay  and  showed  considerable 
protection  at  20  and  2  pM  in  the  MPNHDA.  While  secalonic 
acid  (4),  cephalochromin  (5),  and  skyrin  (6)  were  roughly 
comparable  to  chaetochromin  A  in  the  protease  assay,  they 
were  devoid  of  activity  in  the  MPNHDA. 

In  parallel  with  the  biological  testing,  we  performed  virtual 
testing  of  some  of  these  compounds  using  our  lead 
optimization  tools,^®'^^  including  chaetochromin  A  (l), 
talaroderxine  A  (2),  cephalochromin  (5),  and  skyrin  (6), 
along  with  three  hypothetical  side  chain  oxidized  analogues  of 


Chart  1 

OH  OH  O  OH  OH  O 


1  2 


OH  O  OH 


OH  OH  O 


7H  =  CH2OH,  R2  =  CH3 
8H  R^  =  CH3,  r2  =  CH2OH 
9H  R^  =  R2  =  CH2OH 


chaetochromin  A  (7H— 9H).  Secalonic  acid  was  not  included, 
due  to  substantial  differences  in  its  molecular  scaffold  relative  to 
the  bis-naphthopyrones. 

Because  no  experimental  structural  data  are  available  for 
BoNT/A-bound  states  of  these  compounds,  we  modeled  them 
using  AutoDock  Vina.^^  This  docking  algorithm  preserved 
asymmetric  carbon  chirality  but  was  insensitive  to  the  higher 
order  axial  chirality  of  the  compounds.  Consequently,  our 
modeled  compounds  docked  with  the  BoNT/A  active  site  in 
both  P  and  M  axial  configurations.  The  docking  score  was  not 
sufficient  to  rank-order  closely  related  compounds  and  even 
less  so  to  choose  the  P  atropisomer  over  the  M.  Visual 
inspection  revealed  two  representative  consensus  poses,  one  for 
each  atropisomer,  that  had  carbonyl  and/or  hydroxyl  oxygens 
poised  to  coordinate  active  site  Zn^^.  Figure  1  illustrates  the 
consensus  poses  for  1—3.  However,  only  the  P  atropisomer  was 
capable  of  coordinating  Zn^"^  by  carbonyl  and  hydroxyl  oxygens 
at  the  same  time  (Figure  1A,B).  On  the  basis  of  this  finding  and 
on  the  fact  that  the  original  hit  molecule,  chaetochromin  A  (l), 
had  the  P  configuration,  we  selected  the  pose  with  the  P 
configuration  for  subsequent  binding  free  energy  calculations, 
excluding  talaroderxine  B  (3). 
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Table  1.  Comparative  Testing  of  Fungal  Metabolites“ 


BoNT/A  protease  assay 

%  inhibition 
(/<M) 


mpnhda’’ 

paralytic  time^ 
(/tM) 


sample 

100 

20 

ICjo" 

(/.M) 

20 

2 

chaetochromin  A  (l) 

74.0 

ss.o 

24.6 

187.8^ 

104.0 

talaroderxine  A  (2) 

96.0 

8S.0 

3.7 

197.0^ 

188.7^ 

talaroderxine  B  (3) 

80.3 

67.0 

10.1 

94./ 

174.3^ 

secalonic  acid  (4) 

64.5 

39.1 

28.6 

67.0 

ND^ 

cephalochromin  (S) 

72.0 

4S.0 

29.2 

59.0 

ND^ 

skyrin  (6) 

74.3 

40.3 

31.5 

81.0 

64.0 

CB  7969312^  (control) 

98.0 

96.2 

2.1 

300.0 

300.0 

“Assays  conducted  as  described  in  ref  1  and  the  Supporting 
Information.  *’Mouse  phrenic  nerve  hemidiaphragm  assay.  TC50 
(50%  inhibitory  concentration)  values  were  determined  from  nine 
concentrations  of  each  compound  using  GraphPad  Prism  4  (GraphPad 
Software,  La  Jolla,  CA).  “^Average  time  to  50%  loss  of  twitch  tension 
(min);  average  value  for  BoNT/A  toxin  control,  62  min.  ‘’P  value 
<0.01  (highly  significant)  in  comparison  with  values  recorded  for  the 
BoNT/A  control;  statistical  analyses  were  performed  using  SigmaPlot 
10  (Systat  Software,  San  Jose,  CA).  ^Solubility  issues  in  the  test 
medium  at  this  concentration  rendered  this  value  unreliable.  '^Not 
determined.  ^A  small  molecule  inhibitor  described  in  ref  1. 


Using  the  consensus  pose  with  the  P  configuration,  we  rank- 
ordered  the  ligands  by  their  relative  binding  free  energies, 
which  were  computed  using  computationally  demanding  Single 
Reference  Thermodynamic  Integration  (SRTI)  approach 
augmented  with  Hamiltonian  Replica  Exchange  (HREX).^®’^^ 
Because  binding  free  energy  predictions  often  depend  on  the 
initial  structure  of  the  complex,  it  was  imperative  that  we 
employ  a  similar  pose  for  aU  of  the  compounds  as  the 
consensus  pose.  For  the  same  reason  and  because  axial 
atropisomers  do  not  readily  interconvert  in  SRTI  simulations, 
it  was  not  possible  to  compare  P  and  M  atropisomers  using 
only  the  relative  binding  free  energy  calculations. 

Adequate  sampling  of  the  conformations  of  the  numerous 
hydroxyl  groups  might  be  important  for  obtaining  converged 
results.  Indeed,  each  of  the  studied  compounds  has  at  least  six 
hydroxyl  and  two  carbonyl  groups  that  can  interact  by  forming 
internal  hydrogen  bonds.  The  formation  of  internal  hydrogen 
bonds  would  reduce  the  number  of  favorable  interactions  with 
the  BoNT/A  active  site  and  could  affect  the  binding  free 
energies.  To  ensure  adequate  sampling  of  the  hydroxyl  group 
rotations,  we  employed  a  specially  designed  nonphysical 
reference  state  (see  the  Supporting  Information)  that  allowed 
for  free  rotation  of  the  hydroxyl  groups.  Using  this  reference 
state  in  combination  with  the  HREX  option  within  SRTI 
enhanced  conformational  transitions  of  hydroxyl  groups. 
Therefore,  we  expected  HREX-SRTI  to  provide  more  accurate 
relative  binding  free  energy  predictions  than  conventional  TI 
methods. 

The  results  of  the  calculations  are  summarized  in  Table  2, 
with  the  representative  complexes  of  the  top  scoring 
compounds  shown  in  Figure  1A,B.  Despite  all  of  the 
precautions,  our  predictions  had  relatively  large  standard 
deviations  that,  in  some  cases,  precluded  our  making  a 
definitive  prediction.  Nevertheless,  our  calculations  had 
sufficient  precision  to  predict  talaroderxine  A  (2)  as  a  better 
binder  than  chaetochromin  A  (l).  Interestingly,  one  of  our 
designs  (7H)  was  also  predicted  to  be  a  better  binder  than 
chaetochromin  A.  This  prediction  remains  to  be  tested 


(B)  talaroderxine  A  (2) 


(C)  talaroderxine  B  (3) 


Figure  1.  Proposed  binding  mode  of  (A)  chaetochromin  A  (l),  (B) 
talaroderxine  A  (2),  and  (C)  talaroderxine  B  (3)  in  the  active  site  of 
BoNT/A  The  Zn^"^  ion  is  shown  as  a  magenta  sphere,  carbon  atoms  of 
the  ligands  are  depicted  in  yellow,  and  those  of  the  active  site  residues 
are  in  cyan.  The  molecular  surface  of  the  pocket  is  shown  in  a 
transparent  light  gray  color.  For  clarity,  only  polar  hydrogens  are 
shown.  Panels  A  and  B  depict  the  most  populated  cluster  centers  from 
SRTI-HREX  simulations,  whereas  panel  C  is  the  result  of  docking. 


experimentally,  as  the  compound  is  as  yet  unknown.  Although 
the  precision  of  our  calculations  for  the  remaining  compounds 
did  not  allow  us  to  make  definitive  conclusions,  the  predicted 
changes  in  binding  free  energies  appear  to  be  in  good 
agreement  with  experimentally  determined  IC50  values  (Table 

1). 
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Table  2.  Predicted  Binding  Free  Energies/Standard 
Deviations  Relative  to  Chaetochromin  A“ 


compound 

AAG 

SD 

chaetochromin  A  (l) 

0.0 

l.S 

talaroderxine  A  (2) 

-4.7 

2.1 

cephalochromin  (S) 

1.4 

l.S 

skyrin  (6) 

0.2 

1.1 

7H 

-2.6 

1.3 

8H 

0.7 

1.6 

9H 

-0.3 

2.8 

“Relative  free  energies  were  computed  using  the  Thermodynamic 
Integration  method  and  are  reported  in  kcal/mol.^°’^' 


The  correlation  between  the  predicted  binding  free  energies 
and  the  observed  IC50  values  lends  additional  credibility  to  the 
consensus  pose  selected  for  calculations.  In  the  absence  of  a 
crystallographic  structure,  we  can  use  this  validated  pose  for 
further  lead  optimization  efforts. 

We  should  note,  however,  that  the  higher  activity  of 
talaroderxine  B  relative  to  chaetochromin  A  might  suggest 
that  the  consensus  pose  with  the  M  configuration  is  also 
important  (Figure  1C).  Nevertheless,  the  observation  that 
talaroderxine  A  (2)  is  more  active  than  talaroderxine  B  (3) 
supports  our  emphasis  on  the  P  atropisomer  for  the  present 
study. 

Starting  with  the  initial  in  silico  screening  hit  chaetochromin 
A  (1),  we  compared  a  small  number  of  analogues  and  found 
that  talaroderxines  A  (2)  and  B  (3)  were  more  potent  in  our 
primary  HPLC -based  and  secondary  mouse  phrenic  nerve 
hemidiaphragm  assays.  It  is  noteworthy  that  1  and  2  share  the 
same  axial  configuration  and  absolute  configuration  at  the  chiral 
centers  they  have  in  common;  3  has  the  same  configuration  at 
the  chiral  centers  held  in  common,  but  the  opposite  axial 
chirahty.  The  primary  differences  between  1  and  2  and  3  reside 
in  the  y-  versus  a-pyrone  and  the  methyl  versus  n-propyl 
substituents  at  the  common  chiral  centers.  The  task  before  us 
now  is  to  develop  an  adequate  supply  of  2  and  3,  so  that  we  can 
define  their  mechanism  of  action/mode  of  binding,  assess  their 
in  vivo  activity,  and  determine  whether  they  are  less  toxic  than 
chaetochromin  Because  talaroderxines  A  and  B  (2  and 

3)  are  produced  as  a  mixture  in  relatively  low  yields,  this  will  be 
challenging  but  not  insurmountable.  These  are  the  first  fungal 
natural  products  reported  to  inhibit  BoNT/A 
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